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Closed orbit theory is generalized to the semiclassical cal- 
culation of cross-correlated recurrence functions for atoms in 
external fields. The cross-correlation functions are inverted by 
a high resolution spectral analyzer to obtain the semiclassical 
eigenenergies and transition matrix elements. The method 
is demonstrated for dipole transitions of the hydrogen atom 
in a magnetic field. This is the first semiclassical calculation 
of individual quantum transition strengths from closed orbit 
theory. 

PACS numbers: 32.6Q.-M, 32.70.Cs, 03.65.Sq 



In 1925, Heisenberg ||l| pushed open the door to quan- 
tum mechanics when he proposed a quantum theory 
consisting only of "in principle observable" quantities 
- matrix elements A(n, m) ~, whose physical mean- 
ing was hinged upon the correspondence principle: For 
high quantum numbers n and |n — mj <C n the ma- 
trix elements were to turn into the Fourier amplitudes 
of the corresponding physical observable A{t) associated 
with a classical periodic orbit, viz. A{n, m) e:iqp{i{E{n) — 
E{m))t/h) <^ AT-(n) exp(ira;(n)), with t = n — m and 
a;(n) the frequency of the periodic orbit. By a consistent 
application of the translation rules to the quantization 
condition for actions, AS/ An — h, he arrived at the 
canonical commutation relations, underlying the whole 
of quantum mechanics. Thus establishing the connection 
between the information encoded in classical orbits and 
quantum mechanical transitions amplitudes proved one 
of the fundamental questions of quantum physics. 

The correspondence principle is silent about transition 
amplitudes involving low-lying states. It was, therefore, 
a great success when, more than sixty years later, closed 
orbit theory - a variant of periodic orbit theory H - came 
up with the discovery that there exists an intimate con- 
nection between classical orbits and transitions ampli- 
tudes even in cases where one of the states lies in the deep 
quantum regime. In closed orbit theory, developed by Du 
and Delos |^ and Bogomolny jj], the transition ampli- 
tudes are given as the sum of two terms, one a smoothly 
varying part (as a function energy) and the other a su- 
perposition of sinusoidal modulations. The frequencies, 
amplitudes, and phases of the modulations are directly 
obtained from information contained in the closed classi- 
cal orbits. When the resulting transition amplitudes are 
Fourier transformed the sinusoidal modulations produce 
sharp peaks in the Fourier transform recurrence spectra. 



Closed orbit theory has been applied to the interpretation 
of photoabsorption spectra of atoms in external fields and 
has been most successful in explaining the quantum me- 
chanical recurrence spectra qualitatively and even quan- 
titatively in terms of the closed orbits of the underlying 
classical system [^-0]. 

However, up to now practical applications of closed or- 
bit theory have always been restricted to the semiclassi- 
cal calculation of low resolution spectra for two reasons. 
Firstly, the closed orbit sum requires, in principle, the 
knowledge of all orbits up to infinite length, which are not 
normally available from a numerical closed orbit search, 
and, secondly, the infinite closed orbit sum suffers from 
fundamental convergence problems |g|,y. It is therefore 
usually believed that the calculation of individual tran- 
sition matrix elements, e.g., of the dipole operator D, 
{(j)i\D\ipf), which describe the transition strengths from 
the initial state |0i) to final states \ipf), is a problem 
beyond the applicability of the semiclassical closed orbit 
theory, i.e., is the domain of quantum mechanical meth- 
ods. 

It is the purpose of this Paper to demonstrate that 
high-precision quantum transition amplitudes between 
low-lying and highly excited states can be obtained 
within the framework of closed orbit theory using solely 
the information contained in closed classical orbits. In 
this way we can establish a connection between classi- 
cal orbits and quantum mechanical matrix elements that 
goes far beyond what was conceived of in the early days 
of quantum mechanics. To that end, we slightly general- 
ize closed orbit theory to the semiclassical calculation of 
cross-correlated recurrence functions. We then adopt the 
method of Refs. p|-pJ]|| to harmonically invert the cross- 
correlated recurrence signal and to extract the semiclassi- 
cal eigenenergies and transition matrix elements. Results 
will be presented for the photo excitation of the hydrogen 
atom in a magnetic field. 

The oscillator strength / for the photo excitation of 
atoms in external fields can be written as 



f{E) = (^-^,)Im< 

TT 



\DG+D\(f>,) , 



(1) 



where \ipi) is the initial state at energy Ei, D is the 
dipole operator, and G^ the retarded Green's function 
of the atomic system. The basic steps for the deriva- 
tion of closed orbit theory are to replace the quantum 
mechanical Green's function in (Q) with its semiclassical 
Van Vleck-Gutzwiller approximation and to carry out the 



overlap integrals with the initial state \4'i). Here we go 
one step further by introducing a cross-correlation matrix 



gaa' = {4>a\DG^D\<l)a') 



(2) 



with I0Q,), a — 1,2,...,L a set of independent ini- 
tial states. As will be shown below the use of cross- 
correlation matrices can considerably improve the con- 
vergence properties of the semiclassical procedure. In 
the following we will concentrate on the hydrogen atom 
in a magnetic field (for reviews see [11-13 ) with 7 — 



5/(2.35 X 10'^ T) the magnetic field strength in atomic 
units. The system has a scaling property, i.e., the shape 
of periodic orbits does not depend on the scaling param- 
eter, w = 7~i/^ = h~g , and the classical action scales as 
S = sw, with s the scaled action. As, e.g., in Ref. [|j, 
we consider scaled photoabsorption spectra at constant 
scaled energy E = Ej~^/^ as a function of the scaling pa- 
rameter w. We choose dipole transitions between states 
with magnetic quantum number m = 0. Note that the 
following ideas can be applied in an analogous way to 
atoms in electric fields. Following the derivation of Refs. 
[g|,|| the semiclassical approximation to the fluctuating 
part of gaa' in Eq- 1^ reads 

-(2^)5/2 
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x3^a(^n3^o'(z9}°)e^( 
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the scaled action and Maslov index 
an element of the mon- 



of the closed orbit (co), 
odromy matrix, and t9™ and df the initial and final 
angle of the trajectory with respect to the magnetic 
field axis. The angular functions 3^q(i?) depend on the 
states |(/)a) and the dipole operator D, and are given as 
a linear superposition of Legendre polynomials, 3^0(1?) = 
^; SiQ.P; (cos 1?). For low- lying initial states with prin- 
cipal quantum number n only few coefficients Bia with 
I < n are nonzero. Explicit formulas for the calculation 
of the coefficients can be found in Refs. ||,^. The prob- 
lem now is to extract the semiclassical eigenenergies and 
transition matrix elements from Eq. H because the closed 
orbit sum does not converge. We therefore adopt the 
idea of Ref. ||lj] where we proposed to adjust the Fourier 
transform of a non-convergent Dirichlet series like the 
semiclassical expression (H) to the functional form of its 
quantum mechanical analogue. The Fourier transforma- 
tion of w^^^Qaa' (^) yields the cross-correlated recurrence 
signals 



CO 

with the amplitudes 
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being determined exclusively by closed orbit quanti- 
ties. The corresponding quantum mechanical cross- 
correlated recurrence functions, i.e., the Fourier trans- 
forms of w^''^g^i (w) read 



Caa'(s) = -^^^akha' 



fee 



(6) 



with w-k the eigenvalues of the scaling parameter, and 



bak — 



M^i 



(5) 



"fe yi>a\D\^k) (7) 

proportional to the transition matrix element for the 
transition from the initial state \(^a) to the final state 

The method to adjust (Q) to (||) for fixed states \<\)a) 
and \(i)a') is that of harmonic inversion as discussed in 
Ref. jlj]. However, information theoretical considera- 
tions then yield an estimate for the required signal length, 
Smax ^ 4:Tr g{w) p4| , which may result in an unfavorable 
scaling because of a rapid proliferation of closed orbits 
with increasing period. Moreover, it is a special problem 
to resolve nearly degenerate states and to detect states 
with very low transition strengths from the harmonic in- 
version of a single function C^a'i^)- Note also that the 
element of the monodromy matrix TO12 and the values 
of the angular functions yaii^i) and ya'i'&f) are intrinsi- 
cally intertwined in Eq. |^ for the amplitudes A^a'' i-^-i a 
single function C'aa' (*) does not contain independently 
the information from the monodromy matrix and the 
starting and returning angles di and i?/ of the closed 
orbits. 

In this Paper we therefore propose to apply an ex- 
tension of the method to the harmonic inversion of 
cross-correlation functions [pH10| , which has recently also 
served as a powerful tool for the semiclassical calculation 
of tunneling splittings |15| ] . The idea is that the informa- 
tional content of an L x L time signal is increased roughly 
by a factor of L as compared to a 1 x 1 signal. The ad- 
ditional information is gained with the set of linearly in- 
dependent angular functions yaW, a = 1, 2, ... L in Eq. 
o evaluated at the starting and returning angles i^i and 
■df of the closed orbits. Note that the cross-correlation 
matrix (^) is constructed by using independently the in- 
formation of the closed orbit quantities, i.e., the elements 
77112 of the monodromy matrix and the angles 'di and "df. 
For a given number of closed orbits the accuracy of semi- 
classical spectra can be significantly improved with the 
help of the cross-correlation approach, or, alternatively, 
spectra with similar accuracy can be obtained from a 
closed orbit cross-correlation signal with a significantly 
reduced signal length. 

Here we only give a qualitative and brief description of 
the method. The details of the numerical procedure of 
solving the generalized harmonic inversion problem (|g) 
have been presented in Refs. |@-|l^ . The idea is to recast 
the nonlinear fit problem as a linear algebraic problem [g| . 
This is accomplished by associating the signal Caa' (s) 
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FIG. 1. Real parts of the semiclassical cross-correlated re- 
currence functions for the hydrogen atom in a magnetic field 
at constant scaled energy E = —0.7 as functions of the clas- 
sical action s (in dimensionless scaled units). The graphs of 
the imaginary parts (not shown) qualitatively resemble those 
of the real parts. 

(to be inverted) with a time cross-correlation function 
between an initial state $„ and a final state $„' , 



Cca'{s) = ($a'|e 



— isHf 



«$, 



(8) 



where the fictitious quantum dynamical system is de- 
scribed by an effective Hamiltonian HcS- The latter is 
defined implicitly by relating its spectrum to the set of 
unknown spectral parameters Wk and foafc. Diagonaliza- 
tion of _ffoff would yield the desired Wk and bak- This 
is done by introducing an appropriate basis set in which 
the matrix elements of H^s are available only in terms 
of the known signals Caa'{s). The Hamiltonian HcS is 
assumed to be complex symmetric even in the case of a 
bound system, which makes the harmonic inversion sta- 
ble with respect to "noise" due to the imperfections of 
the semiclassical approximation. 

We now demonstrate the method of harmonic inversion 
of the cross-correlated closed orbit recurrence functions 
(0) for the example of the hydrogen atom in a magnetic 
field at constant scaled energy E = —0.7. This energy 
was also chosen for detailed experimental investigations 
of the helium atom 0] . We investigate dipole transitions 
from the initial state \(j)i) = |2p0) with light polarized 
parallel to the magnetic field axis to final states with 
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FIG. 2. Quantum (upper part) and semiclassical (lower 
part) photoabsorption spectra of the hydrogen atom in a mag- 
netic field at scaled energy E — —0.7 as functions of the di- 
mensionless scaling parameter w — 7"^'"': Transition matrix 
elements {2pO\D\ipf)^ for dipole transitions with light polar- 
ized parallel to the magnetic field axis. 

magnetic quantum number m = 0. For this tran- 
sition the angular function in Eq. [^ reads 3^1 (i?) — 
(27r)-i/22'^e"''(4cos2i? - 1). For the construction of a 
2x2 cross-correlated recurrence signal we use for simplic- 
ity as a second transition formally an outgoing s-wave, 
i.e., D\(j)2) oc Yofi, and, thus, 3^2(1?) — const. A nu- 
merical closed orbit search yields 1395 primitive closed 
orbits (2397 orbits including repetitions) with scaled ac- 
tion s/2tt < 100. With the closed orbit quantities at hand 
it is straightforward to calculate the cross-correlated re- 
currence functions in (H) . The real parts of the functions 
Cll{s), 012(3), and C|2(s) with s/27r < 50 are presented 
in Fig. n^. The imaginary parts are not shown because 
they qualitatively resemble the real parts. Note that for 
symmetry reasons 6*21(5) — Cf2(s). 

We have inverted the 2x2 cross-correlated recur- 
rence functions in the region < s/2tt < 100. The 
resulting semiclassical photoabsorption spectrum is com- 
pared with the exact quantum spectrum in Fig. ^a for 
the region 16 < w < 21 and in Fig. ^b for the region 
34 < w < 40. The upper and lower parts in Fig. g 
show the exact quantum spectrum and the semiclassi- 
cal spectrum, respectively. Note that the region of the 
spectrum presented in Fig. j^b belongs well to the experi- 



mentally accessible regime with laboratory field strengths 
B = 6.0T to B = 3.7T. The overall agreement between 
the quantum and semiclassical spectrum is impressive, 
even though a line by line comparison still reveals small 
differences for a few matrix elements. It is important 
to note that the high quality of the semiclassical spec- 
trum could only be achieved by our application of the 
cross-correlation approach. For example, the two nearly 
degenerate states at w; = 36.969 and w = 36.982 cannot 
be resolved and the very weak transition at w = 38.894 
with (2pO|I?|V'/)^ = 0.028 is not detected with a sin- 
gle (1x1) recurrence signal of the same length. How- 
ever, these hardly visible details are indeed present in the 
semiclassical spectrum in Fig. |2fD obtained from the har- 
monic inversion of the 2x2 cross-correlated recurrence 
functions. 

In conclusion, we have demonstrated that closed orbit 
theory is not restricted to describe long-range modula- 
tions in quantum mechanical photoabsorption spectra of 
atoms in external fields but can well be applied to extract 
individual eigenenergies and transition matrix elements 
from the closed orbit quantities. This is achieved by a 
high resolution spectral analysis (harmonic inversion) of 
cross-correlated closed orbit recurrence signals. For the 
hydrogen atom in a magnetic field we have obtained, for 
the first time, transition matrix elements between low- 
lying and highly excited Rydberg states using exclusively 
classical closed orbit data. It will be straightforward, and 
rewarding, to apply the method to atoms in electric fields. 
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